Diffusion Limited Aggregation and Iterated Conformal Maps 
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The creation of fractal clusters by diffusion limited aggregation (DLA) is studied by using iterated 
stochastic conformal maps following the method proposed recently by Hastings and Levitov. The 
object of interest is the function $^ n ' which conformally maps the exterior of the unit circle to the 
exterior of an n-particle DLA. The map $' n ' is obtained from n stochastic iterations of a function 
4> that maps the unit circle to the unit circle with a bump. The scaling properties usually studied 
in the literature on DLA appear in a new light using this language. The dimension of the cluster is 
determined by the linear coefficient in the Laurent expansion of & n > , which asymptotically becomes 
a deterministic function of n. We find new relationships between the generalized dimensions of the 
harmonic measure and the scaling behavior of the Laurent coefficients. 



I. INTRODUCTION 

The diffusion limited aggregation (DLA) model was in- 
troduced in 1981 by T. Witten and L. Sander §J. The 
model has been shown to underlie many pattern form- 
ing processes including dielectric breakdown M , two-fluid 
flow pj , and electro-chemical deposition Q . The model 
begins with fixing one particle at the center of coordi- 
nates in ci-dimensions, and follows the creation of a clus- 
ter by releasing random walkers from infinity, allowing 
them to walk around until they hit any particle belong- 
ing to the cluster. Upon hitting they are attached to 
the growing cluster. The model was studied on- and off- 
lattice in several dimensions d > 2; here we are only in- 
terested in the off-lattice versions in two dimensions. 

DLA has attracted enormous interest over the years 
since it is a remarkable example of the spontaneous cre- 
ation of fractal objects. It is believed that asymptotically 
(when the number of particles n — » oo) the dimension D 
of the cluster is very close to 1.71 ||, although there ex- 
ists to date no accepted proof for this fact in spite of 
several interesting attempts [[^j . In addition, the model 
has attracted interest since it was among the first || to 
offer a true multifractal measure: the harmonic measure 
(which determines the probability that a random walker 
from infinity will hit a point at the boundary) exhibits 
singularities that are usefully described using the multi- 
fractal formalism ||. Nevertheless DLA still poses more 
unsolved problems than answers. It is obvious that a new 
language is needed in order to allow fresh attempts to ex- 
plain the growth patterns, the fractal dimension, and the 
multifractal properties of the harmonic measure. 

Such a new language was proposed recently by Hast- 
ings and Levitov Q. |ilj| . These authors showed 
that DLA in two dimensions can be grown by iterating 
stochastic conformal maps. We adopt their basic strat- 
egy and will see that it provides a new formulation of the 
problem which may lead to new insights and results. 



The basic idea is to follow the evolution of the confor- 
mal mapping (w) which maps the exterior of the unit 
circle in the mathematical w-plane onto the complement 
of the cluster of n particles in the physical z-plane. $W 
is unique by the Riemann mapping theorem, provided 
that it satisfies the boundary condition 



r?( n ) 



as w 



(1) 



Here F[ n) is a real positive coefficient, fixing the argu- 
ment of [<&( n >(w)]' to be zero at infinity. <S>( n \w) is re- 
lated to the complex electric potential ^ n \z) by 



¥ n \z) = lnM")(z) : 



(2) 



where hS n '{z) — (z) is the inverse mapping. Let- 

ting z — > oo in Eq.(|l|) it is easy to verify that Eq.(||) 
implies 



^ (n \z) — lnz when 



oo 



(3) 



as it should be at d = 2. 

The equation of motion for $' n J (w) is determined re- 
cursively. The choice of the initial map $ < -°- ) (w) is rather 
flexible, and in this paper we select (arbitrarily) an ini- 
tial condition <E>(°)(ui) = w. We expect the asymptotic 
cluster to be independent of this choice. Then suppose 
that <I>( n ~ 1 ^(w) is given. The cluster of n "particles" is 
created by adding a new "particle" of constant shape and 
linear scale a/Ao t o the cluster of (n — 1) "particles" at a 
position which is chosen randomly according to the har- 
monic measure. We denote points on the boundary of the 
cluster by Z(s) where s is an arc-length parametrization. 
The probability to add a particle on an infinitesimal arc 
ds centered at the point z(s) on the cluster boundary is 



P(s,ds) ~ |W(s)|ds. 



(4) 



The pre-image of z(s) and ds in the w-plane are e and 
dO respectively. Clearly, ds = \ [& n -V}'(e ie )\d6. From 
Eq.(0) we conclude that 



1 



P{s,ds) ~ |W(s)||$'|d0 = d6, 



(5) 



so the harmonic measure on the real cluster translates to 
a uniform measure on the unit circle in the mathematical 
plane. 

The image of the cluster of n particles under (z) 
is, by definition, just the unit circle. On the other hand, 
the image of the cluster of n particles under h!- n ~ x \z) 
is the unit circle with a small bump whose linear scale 
is \/Ao/|^™ (e lSn )\ where e l9n is the image (under 
of the point z n on the real cluster at which the 
growth occurred. 

Let us define now a new function 4>\ n .g n (w). This func- 
tion maps the unit circle to the unit circle with a bump 
of linear scale y/X^ around the point e lBn . For w — ► oo, 
(j)\ n e n {w) ~ w (with positive real proportionality coeffi- 
cient). Using </>A„,e„ (u>) the recursion relation for $l n J (w) 
is given by (see Fig. |l|): 



According to the above discussion A n is given by 

A 



An 



|$(n-l)'( e ie„)|i 



(6) 



(7) 



so the RHS of Eq.(^|) is determined completely by 
(w); Eq.@ induces the recursive dynamics of 
$(")(«;). 

The recursive dynamics can be represented as itera- 
tions of the map (f>x„,e n (w), 

$ (n) (w) = (f)\ 1 i g 1 o <j) X2fi2 o . . . o «/>A„,e„M ■ (8) 

This composition appears as a standard iteration of 
stochastic maps. This is not so. The order of iterations is 
inverted - the last point of the trajectory is the inner ar- 
gument in this iteration. As a result the transition from 
<f>( n \w) to $( n+1 )(u;) is not achieved by one additional 
iteration, but by composing the n former maps Eq.(|^) 
starting from a different seed which is no longer u> but 

0A„ +1 ,0„ +1 (w)- 

We note that in the physical plane the "particles" are 
roughly of the same size. To achieve this the linear scales 
VA^ vary widely as a function of n and 9. We will see that 
the distribution of y^n and their correlations for differ- 
ent values of n determine many of the scaling properties 
of the resulting cluster. In particular their moments are 
related to the generalized dimensions of the harmonic 
measure. 

There are many functions <p\,e which conformally map 
the unit circle to the unit circle with a bump. A simple 
choice is a function which behaves linearly for large w and 
has a simple pole inside the unit circle which will induce 
a bump in the image. The pole has to be at wq = 1 — A 
in order to localize the bump near w = 1 and make it of 
linear size of the order VA. The residue has to be A 3 / 2 , 



in order for the bump's height to be also of the order 
Consider then 



(w) = (1 + X)w + 



A 3/2 



w — wo 

Careful thinking leads to the conclusion that this func- 
tion and other similar functions are inappropriate: they 
have long "tails" . In other words, the unit circle is 
slightly distorted everywhere. This small global distor- 
tion may result in a loss of conformality or in the growth 
of non-constant size particles in the physical plane in nu- 
merical applications. 

It was proposed in Ref. ]l0| that a choice for 4>x n ,e n (w) 
that is free of global distortion is given by 



1 + w + w[ 1 + 



2 iv 
1 



2 1 - A 
wl + A 



1/2' 



- 1 



<Px,e(w) = e ie (f>xfi(e~ ie w) , 



(9) 
(10) 



The parameter a is confined in the range < a < 1. As 
a decreases the bump becomes flatter, with the identity 
map obtained for a — 0. As a increases towards unity 
the bump becomes elongated normally to the unit circle, 
with a limit of becoming a line ( "strike" in the language 
of jl0|]) when a = 1. Naively one might think that the 
shape of the individual particle is irrelevant for the large 
scale fractal statistics; we will see that this is not the 
case. The dependence on a is important and needs to 
be taken into account. Notice that this map has two 
branch points on the unit circle. The advantage of this is 
that the bump is strongly localized. On the other hand 
repeated iterations of this map leads to rather complex 
analytic structure. 

The aim of this paper is therefore to investigate the 
scaling and statistical properties of such iterated stochas- 
tic conformal maps with a view to discovering the scaling 
properties induced by the dynamics which any analytic 
theory must ultimately explain. In Section II we present 
the numerical procedure used to generate the fractal clus- 
ters, and in Section III give the necessary mathematical 
background to describe such mappings. In particular we 
discuss the Laurent expansion of the conformal map from 
the unit circle to the n-particle cluster; the coefficients 
of the Laurent series have interesting scaling behavior 
with the size of the cluster which is intimately related to 
the fractal dimension of the cluster and to the general- 
ized dimensions of the harmonic measure. In Section IV 
we present numerical results regarding the scaling prop- 
erties of averages of the Laurent coefficients and of the 
size parameter A„. The results are accompanied by a 
theoretical analysis and interpretation. In Section V we 
conclude with some remarks on the road ahead. 
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II. NUMERICAL PROCEDURE 

The algorithm simulating the growth of the clus- 
ter is based on Ref. |ic}| . The n "particle" cluster 
is encoded by the series of pairs {(8i, Ai)}™ =1 . Hav- 
ing the first n — 1 pairs, the n th pair is found as 
follows: choose 9 n from a uniform distribution in 
[0, 27r], independent of previous history. Then com- 
pute A„ from Eq.(^j), where the derivative of the it- 
erated function involves 4>'\ _ g _ 1 , 4>'\ 7l _ 2 e„_ 2 > 
4>'x _ e e ^ c ' computed respectively at the points 

e»ViL S ll «»-i(e W ") J ^-■.fln-^An-Lfln-^e*")), etc. 
Notice that the evaluation of both (j)' and cj> after the 
addition of one particle involves 0(n) operations since 
the seed changes at every n. This translates into n 2 time 
complexity for the growth of an n-particle cluster. This 
is inferior to the best algorithms to grow DLA (using hi- 
erarchical maps , with close to linear efficiency) , but 
the present algorithm is not aimed at efficiency. Rather, 
it is used since the Laplacian field and the growth prob- 
ability which is derived from it are readily available at 
every point of the cluster and away from it. The typical 
time to grow a 10,000 particle cluster is 8 minutes on a 
300 MHz Pentium-II. 

Naively one would expect that any choice of < a < 1 
would yield DLA clusters, since a only determines the 
shape of the particles (the aspect ratio is \a/{l — a) for 
small A), and the microscopic details of the particles (ex- 
cept their linear size) should not affect the global prop- 
erties. Three typical clusters with particles of various as- 
pect ratios a are shown in Fig. We mark in black the in- 
terior of the image of the unit circle under the conformal 
map & n \w). The objects look very much like typical 
DLA clusters grown by standard off-lattice techniques, 
and in the next section we demonstrate that they have 
fractal dimensions in close agreement with the latter. For 
a significantly different from 2/3, disadvantages of the al- 
gorithm get amplified. Since the functional form of 4> is 
fixed (only the size and position of the "bump" change) , 
particles of constant shape and size are obtained only 
if the magnification factor |$( ra_1 ) | (the inverse of the 
field) is approximately constant in the iw-plane around 
the "bump" of (f>. If the particles are elongated along 
the cluster, then the variation of the field along the clus- 
ter affects the shape: large otherwise deeply invaginated 
regions, where <!>' is large, are filled up with a single par- 
ticle, and the resulting cluster tends to be more compact. 
This effect, slightly noticeable even at a = 2/3, is quite 
significant at the otherwise natural choice of a = 1/2, 
where the particles are half circles. In Fig. || we show 
such a cluster and point out to the area filling dark re- 
gions which represent such unwanted events. The other 
extreme, when the particles stick out of the cluster, leads 
to sensitivity to variations in the field going away from 
the cluster. Especially if a bump is grown on a tip of 
a branch, where the field decreases rapidly as one goes 
away from the tip (such that increases significantly), 



then the map of the bump gets magnified, resulting in 
particles of very unequal sizes. 

It is necessary to stress that even for a — 2/3, when 
this procedure appears to yield nice ramified structures, 
the problem of fill-ups does not go away: in a few rare 
cases the particle - if it happens to land on a place where 
|<fr"| is large - is significantly distorted. The net effect is 
that large areas surrounded by the cluster (where the 
growth probability is small) are filled up entirely by one 
distorted particle. For the value of a = 1/2 it appeared 
that the errors may be unbounded. Our numerics in- 
dicates that for a = 2/3 the errors were bounded for 
the cluster sizes that we considered. We do not have 
a mathematical proof of boundedness of the errors, and 
our disregard of this danger is only based on the sensible 
appearance of our clusters at this value of a. 



III. MATHEMATICAL BACKGROUND 

In this section we discuss the Laurent expansion of 
our conformal maps, and introduce the statistical objects 
that are studied numerically in the next section. 



A. Laurent Expansion 

Since the functions $>( n \w) and cj>\^(w) are required 
to be linear in w at infinity, they can be expanded in a 
Laurent series in which the highest power is w: 

$(«) ( w )= F^w + F (n) + F^w- 1 + F { ^w- 2 + ... (11) 
M™)= hw + fo + f-iw- 1 + f- 2 w- 2 + ... (12) 

where 

fi = (1 + A) Q 

, 2a\e ie 
Jo - 



f-i 



(l + A) 1 -" 

2aXe 2W 
(1 + A) 2 -" 

a 3i9 



2a- 1 , 
1 H — A 



2aXe m ( 

<£M 1+a( *- 1,A+ 



2a 2 - 3a + 1 2 



The recursion equations for the Laurent coefficients of 
can be obtained by substituting the series of <3? 
and d> into the recursion formula (p|). We find 



F (n) 


= F[ n ~ 


-1) An) 
Jl 


F (n) 
r 


= F[ n ~ 


l)f(n) 
■>0 


F (n) 


= F[ n ' 


-1) j,(n) 


F (n) 


= i>~ 


-1) An) 
■1-2 



? (n-l) 


? (n-l) i „(n) 



f{n) 



F ( 



(/ (") )2 



F 



(n-l) 



(/ (n) )2 



(13) 



(14) 
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We note that the n-dependence of follows from the 
dependence on the randomly chosen 9 n at the nth step, 
from which follows the dependence of A n on n. The lat- 
ter is however a function of all the previous growth steps, 
making the iteration ( |l3| ) -((Til) rather difficult to analyze. 

A general relation between the Laurent coefficients is 
furnished by the so-called area theorem which applies to 
univalent mappings. Since our maps solve the Laplace 
equations with boundary conditions only at infinity and 
on the cluster boundary where the potential is zero, they 
map the w plane uniquely (and with a unique inverse) to 
the z plane. In other words, the pressure lines and the 
stream lines are non-degenerate. Such mappings have the 
property [jl3| that the area of the image of the unit disc 
in the nth step is given by: 



F 



O) 



k=l 



(15) 



A second theorem that will be useful in our thinking is 
a consequence of the so-called one-fourth theorem, see 
Appendix A. There a statement is proven that the inte- 
rior of the curve {z : z = <j>( n ) (e 10 )} is contained in the 

z-plane by a circle of radius AF[ n ^ . Now as the area S n 
is obtained simply from the superposition of n bumps of 
roughly the same area Aq, it has to scale like S n » nA , 
for large n. On the other hand any typical radius of the 
cluster should scale like n^^^/Ao where D is the dimcn- 

(n) 

sion of the cluster. We can thus expect a scaling of F± 
that goes like 



(n) 



(16) 



We note in passing that this scaling law offers us a very 
convenient way to measure the fractal dimension of the 

growing cluster. Indeed, we measured the dimension D 

(n) 

for a range of a in this way by averaging F± over 100 
clusters. We found that for a range of a spanning the 
interval [1/3,8/9] the dimension is constant, around 1.7. 
We can infer from Eq.(fl6|) that the sum in Eq. jl5| ) 

which subtracts positive contributions from |F 1 ^™' ) | 2 con- 
tains terms that cancel the behavior of n 2 ^ D (remember 
that D < 2), leaving a power of unity for the scaling of 
S n . Indeed, we will show below both numerical and the- 
oretical evidence for the scaling behavior of the l-fl^l 2 

for k > 6 which is in agreement with n 2 ' D . 

We can give a direct physical interpretation for the co- 
efficients F^' by comparing them to the coefficients of 
the series for cf. Eq.(§): 



i z 



(17) 



The coefficient of ln(z) is unity so that the electric flux is 
unity. This corresponds to the normalization of the prob- 
ability. The constant ro is the Laplace radius which is the 



radius of a charged disk which would give the same field 
far away. The rest of the V'fc's are conventional multipole 
moments. 

The relations between the Laurent coefficients of \Jj( n ) 
and are: 

r = F t 



lfj 4 



= -Fa 



-F^Fi - -Fl 

-F-2F 2 — 2F F-iFi — -Fq 

-F_ 3 F? - -F^Ft - 3-Fi-F 2 -F_i 



-3F- a F F} - -F* 



(18) 



The first line shows that F% = ro, the Laplace radius, 
in accordance with the one-fourth theorem . 

The second line shows that the dipole moment ipi is 
— Fq. We can interpret this coefficient as a distance, the 
wandering of the center of charge due to the random ad- 
dition of the particles. We will take the point of view that 
this quantity is less "intrinsic" than the others to the dy- 
namics of the DLA growth. In fact, if we set Fo = ipi = 0, 
(we could imagine shifting the cluster as we grow it) we 
can rewrite the rest of the equations: 



-F-a 



(V>4 + \^D/<\ 
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• 



(19) 



etc. This leads to the interpretation of in terms of 
the multipole expansion of the electric field. 



B. Statistical objects and the relations to 
generalized dimensions 

Our growth process is stochastic. Accordingly, it is 
natural to introduce averages over the randomness. In 
our thinking there are two important averages, one over 
histories of the whole random trajectory {0i}f_ 1: and the 
other only over the random choice of 6 n at the nth step. 
To distinguish between the two we denote the first by an- 
gular brackets and refer to it as "history-average" , while 
the second is denoted by an overbar and referred to as 
a "cluster-average" . There is a possibility that for very 
large clusters (n — * 00) the two averages result in the 
same numbers. We will refer to such a property as "self- 
averaging" . 

The cluster average of moments of A„ offers a relation- 
ship to the generalized dimensions of the harmonic mea- 
sure [Q. The latter are defined by dividing the plane 
into boxes of size e, and estimating the probability for a 
random walker to hit the piece of the boundary of the 
cluster which is included in the zth box by 



4 



Pile) 



\Ei\e 



(20) 



A. Laurent Coefficients of $ (n) 



where \Ei\ is the modulus of the electric field at 
some point in the ith box. The generalized dimensions 
arc defined by the relation 



N(e) 



(21) 



where N(e) is the number of boxes of size e that are 
needed to cover the boundary, and R is the linear size 
of the largest possible box, which is of the order of the 
radius of the cluster. Substituting d2l]) we find 



N(e) 
1=1 



(s-i)A, 



(22) 



Taking e very small, of the order of \/Ao, and assuming 
that the field is smooth on this scale we have: 



\E^ds~ (VA 



,1-9 n Q-<l)DJD 



(23) 



where L is the length of the boundary, ds is an arc- 
length differential, and we have used the scaling law 

The connection to our language is obtained by con- 
sidering the cluster average of powers of A„. We grow 
a cluster of n — 1 particles, perform repeated random 
choices of growth sites (without growing) , and compute 
X n for each choice. The cluster average can be repre- 
sented as an integral over the unit circle, A^, and is given 
by 

XI = (1/2tt) / Xl(d)dd . (24) 
Jo 

Recalling Eq. (@) we observe that A„(0) = Xl\E(6)\ 2q . 
The last relation, Eq.(||), and Eq.(|2^) imply the scaling 
relation 



An 



-2qD 2q+1 /D 



(25) 



IV. 



NUMERICAL RESULTS AND THEIR 
INTERPRETATION 



In this section we present results on three topics: 

(i) The coefficients of the Laurent expansion. The scaling 
behavior of these quantities is described and discussed in 
the first subsection. 

(ii) The microscopic fluctuations in the conformal map. 
We show that the assumption of self-averaging is valid 
for Eq. ( p5| ) and that the multi-fractal exponents are in 
a good agreement with the known ones. 

(iii) Distribution functions of the Laurent coefficients. 
We analyze numerically the width of those functions and 
find that tends to a deterministic function of n. We 
attribute this effect to non-trivial temporal correlations 
in the field, and give some evidence of their existence. 



All the coefficients of the Laurent series of $(»)(ty) are 
complex numbers except F\ which is real by the choice 
of zero phase at infinity, see Eq. ([!]). Most of our discus- 
sion below pertains to the amplitudes of the coefficients 
Ffe. We need to stress, however, that the phases are not 
irrelevant. If we attempted to use the correct amplitudes 
with random phases, the resulting series will in general 
not be conformal. 

One of the main results of this paper is that in a addi- 
tion to the expected scaling behavior of the linear coeffi- 
cient F[ n) ( given in Eq. (|l6|)) the rest of the amplitudes 
of the Laurent coefficients l-F^'l exhibit also a scaling 
behavior. We find numerically that in the mean the mag- 
nitudes of the Laurent coefficients scale as powers of n: 



a k n 



(26) 



The exponents x k are given in Fig. a. We first discuss 
the consequences of the scaling behavior of F± . 



1. Scaling of Fi 

The scaling behavior (|l6| ) has immediate consequences 
for the scaling behavior of the bump areas A n . The con- 



nection appears from the recursion Eq. ( |13| ) of F^ 
together with /i = (1 + A) a reads 



(n) 



F^ = l[[l + X k } a . 



which 



(27) 



k=l 

Taking history averages we find 

n 

k=l 

n 

ln(F 1 ( " ) )«a^(A fc ) 
fe=i 



d\n{F[ n) )/dn^a{X n ). 



(28) 

(29) 
(30) 



The last two equations are obtained by expanding the 
logarithm and keeping only divergent sums: both the 
mean of F^ and the mean of the sum of Afc increase as 
a function of n. All other sums of powers of Afc converge 
as a function of n, cf. subsection B. Thus, if we assume 
that (F^ n) ) (x n 1 / , cf. Eq.([j]|K fractal scaling of the 
radius (see below) implies that |10(] 



(A n ) = 1/naD. 



(31) 
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In the next subsection we show that this is indeed sup- 
ported by the simulations. Note that (A„) is inversely 
proportional to n for any value of the fractal dimension 
D. On the other hand, if we assume the property of self- 
averaging, Eq. implies a relationship between the 
generalized dimension D3 and the fractal dimension D. 
Comparing Eqs.(p5j) and ( pl| ) leads immediately to the 
relation 



D 3 = D/2 . 



(32) 



This scaling relation was derived by Halsey |l5| using 
much more elaborate considerations. We see that in the 
present formalism this scaling relation is obtained very 
naturally. In fact the present formulation is more pow- 
erful since Eq. ( |3l| ) predicts not only the exponent of the 
third moment of the electric field, but also the prefactor. 
It is also noteworthy that the scaling relation (J3^) results 
simply from the existence of a power law behavior for the 
radius F± 



(<>) 



2. Scaling of Fq 

We found the exponent of (|i*b| 2 ) to be x — 0.7 ± 0.1, 
see Fig. f|. To estimate the scaling behavior of F theo- 
retically we note that 

Fo = * {n){e)d6 = ^ £ z ( s )\ E ( s )\ ds ■ (33) 



Accordingly we can write 

F, -' = (1/47T 2 ) 



- - ' i ^ / ds / ds'z(s)z{s')*\E(s)\\E{s / )\ 
Jo Jo 

L 

2 / j„Ie>/„m2 



X a R z I ds\E{s)\* 
lo 



(34) 



In writing the second line we assumed that the main con- 
tribution to the correlation function is short ranged, 



(z(s)z(s'y\E(s)\\E(s')\) ~ XoR 2 \E( s )\^( t 



(35) 



The justification for this is that the field is expected to 
exhibit wild variations as we trace the boundary z(s). In 
addition the main contribution to the integral is expected 
to come from the support of the harmonic measure where 
the radius is of the order of R. ^From the estimate (|34|) 
and Eq.© we then find 



2- D 2 
x = » 0.64 



(36) 



in agreement with our measurement of xq ■ (We used here 
Z?2 = 0.90 in correspondence with the numerical finding 
reported in Section IV C| . Any of the values of Di quoted 
in the literature would yield xo in the range 0.7 ± 0.1.) 



3. Scaling of F_ k 

The exponents Xk for k < are smaller than 2/D but 
approach it asymptotically, see Fig. [| This behavior is 
expected from the area theorem, and also from a direct 
estimate of the integral representation of the coefficients 
for large k 



1 



-fc = 



ds 



x / ds'z(s)z( S ')*\E(s)\\E( S ')\e^ e ^- e ^''> . (37) 

In Appendix B we show that this integral can be esti- 
mated using the multifractal formalism of the harmonic 
measure with the final result 



{R/Akf / da(2k/Tr) f{a)/a 



(38) 



where a and f(a) are the strength of singularities of the 
harmonic measure and the dimension of the sets of points 
that exhibit these singularities respectively |BL For our 
purposes the important consequence of Eq.(38l) is the 
scaling relation (assuming self- averaging) 



(\F_ k \ 2 )=\ Q n 2 / D g(k) 



(39) 



with g(k) - 1/k 2 J dak^ a ^ a . One knows fr om the 
theory of multifractals that f(a)/a < 1, and there- 
fore we can bound g(k) from above and from below, 
Akr 2 < g(k) < Bk~ x . This is in accord with our numer- 
ical simulations in the range 3 < k < 10, although the 
calculation in the appendix is only valid for large values 
of k. We found agreement with Eq.(p6|) with x k — > 2/D 
and a/c ~ k~ a with 1 < a < 2 . 

Note that this scaling behavior has important conse- 
quences for both the area theorem and for conformality. 

Absolute convergence of the sum J^kLi H^-k I 2 m tne 
area theorem requires a > 2 which is not the case. The 
situation is even more serious for the existence of con- 
formality. To insure the latter the sum J2kLi M^-k I 
must exist. This would require a > 4. The reason 
that the sums exist in the theory is only due to the ul- 
traviolet cutoff at vAo- This cutoff introduces a high- 
est k in the Laurent expansion which we estimate as 
2Trk max i=s L/yAo ~ ti where L is the perimeter of the 
cluster. 



B. Multi- fractal exponents 

Here we test Eq.(p5|). In Fig. || we display double- 
logarithmic plots of (A£) vs.n for q — 0.5, 1, 1.5, 2, 2.5, 3 
and 3.5. The values of the exponents derived from our 
simulations agree very well (within the uncertainties) 
with the generalized dimensions D q obtained in the past 
H for D 2 , ■ ■ ■ , D s using standard methods. In addition 
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we reproduce numbers in agreement with the theoretical 
prediction of D = D « 1.71 and D 3 = D/2. This 
agreement is a strong indication for self averaging at 
least for the purpose of computing moments of A„ (i.e. 
<A<0~Al). 



C. Fluctuations of the averages 

We previously discussed the scaling behavior of \F^ 
and showed that their history averages obey Eq. ( 26[ ) 
However \F 



'_ k | are random variables with broad scaling 
distributions. Fig. |^ describes the rescaled standard de- 
viation cr[™' ) of the Laurent coefficients, 



» 



IF, 



(n)i 



/(IF; 



(n)i 



(40) 



for k — 1,0,-1,-2 as a function of the cluster size n. 
As is seen clearly from the graphs the widths of the dis- 
tributions for all k < tend asymptotically to a finite 
value. This is the normal behavior for scaling distribu- 
tions. The exceptional case is k = 1. Even though it 
exhibits a scaling law of the type (^) (see Section III), 
with 



. Tl = _~1.18, 



the rescaled distribution width of \F^\ tends to zero 
as n goes to infinity. This means that the rescaled distri- 
bution function of tends asymptotically to a delta 
function. The importance of this result for the evalua- 
tion of the fractal dimension of the cluster warrants an 
immediate discussion of this sharpening phenomenon. 

The conclusion of the numerics on F\ is that there ex- 
ists a universal constant c(Aq) such that 



-1/Dp(n) 



c(A ) 



(41) 



where c(Ao) is cluster independent! Moreover, we found 
that c(Ao) = c\f\$, which is in accordance with the role 
played by as an ultraviolet inner lengthscale, which 
is the only lengthscale that appears in the mappings. 
Note that the constant c in Eq.(|4l]) depends on the pa- 
rameter a. We measured c values of 0.6, 0.87, 1.2 and 1.8 
for a values of 1/3, 1/2, 2/3 and 4/5 respectively. 

The observed sharpening is not obvious since we 

(n) 

knowthat F± is built from a product of random vari- 
ables A„, whose moments change with n in multi- fractal 
manner according to Eq. (|25|). 

One could attempt to connect the sharpening of F± 
to the existence of other sharp functions of n. Consider- 
ing the full expansion of Eq. (E7h we find 



1 



(42) 



n 1 n 1 n 

e^-Je^e^ 



2^ 1 3 

i=l i=l 



(43) 



We could understand Eq. (^l|) easily if all the sums of all 
the powers of \ converged to constants, 



VA 4 --ln 
D 



n — > ci 



E A * 



C-2 



i = l 



(44) 

(45) 
(46) 



with Ci cluster independent. In fact, this is not the case. 
The sums of powers are not cluster independent. A clear 
demonstration of this is a simulation that we performed 
in which the initial condition was very far from a circle. 
The individual sums in Eq. ( 43 ) were very different from 
the average values, but nevertheless m (l + con ~ 
verged to the right value. It is our conclusion that each 
of the sums in ([l3|) is not cluster independent, and yet 
somehow the resummed form is cluster independent. 

This remarkable sharpening calls for further discus- 
sion; it appears that its interpretation requires better 
understanding of the time correlations of the field: an 
independent choice of random realization of a series of 
Xi according to their multi-fractal distribution can only 
generate F^ with the proper scaling exponent but can- 

(n) 

not trivially yield a highly peaked distribution of F± . 
Therefore we consider now some evidence for the exis- 
tence of temporal correlations. 

The first outstanding evidence appears in the context 
of the scaling behavior of Fo, which was discussed in the 
first subsection. We show that if we assume that there 
exist no correlations between different growth stages, the 
exponent xq will be very different from the measured 
and calculated value. ^From the recursion relations of 
the Laurent coefficients (Eq. |l4|) we can estimate, in the 
limit of large n when A n is very small on the average, 



|F n ^ 



£ E (^ (m) ^ (m ' 3 A m A r , 



(47) 



m=l m' — l 
n n 



J- J2 (F^F^)(\ m \ m ,eW™- e ^) (48) 



The second line is obtained because F^ is proportional 
to the radius of the whole cluster and should not be cor- 
related with A m . The crucial approximation comes next: 
if A m and X rn i can be treated as independent for m =/= m' , 
then (since 9 m and # m < are independent) Eq.(^8|) simpli- 
fies to 



(\Ft ] \ 2 )- 

n 



(49) 



(50) 
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The numerical simulation resulted in an exponent of the 
order of 0.7, in serious disagreement with Eq. ([5Ch. We 
think that the assumption of independence, Eq. (j49|) is 
the culprit. 

Another fact which illustrates the importance of the 
time-angle correlation Eq. ([l9]) is the difference between 
the exponents of Fo and F-i ((|Fo| 2 ) ~ n - 7 whereas 
(|FLi| 2 ) ~ n 9 ). Their equations of motion ( |l4]) dif- 
fer, for small A n , by two terms only. The first one is 
the term A„F^" _1 ' ) in the RHS of the equation for F_i 
which is absent in the equation for Fo. We checked nu- 
merically that neglecting this term leads to a very small 
change in the exponent. The second difference between 
is that the term A n A,] 



_ fee H0»-0»-*J i n Eq.(^8|) is replaced 
by A„A„_fce 2l( - t '™ _w ' l -'=). The change in the exponent can 
therefore be directly attributed to the existence of im- 
portant time-angle correlations. 

We tried to analyze numerically the time-angle corre- 
lations (A n A„_ / fce l ( e ™~ e '™- fc )). The results for some fc's are 
shown in Fig. ^. It appears that as we increase the size 
of the ensemble, (A n A„_fee l ^ ,, ~ e '*- fc )) — > with the usual 
N^ 1 / 2 dependence on the ensemble size. If we believe 
these numerical results (doubts may exist due to the rel- 
ative smallness of the ensembles analyzed), then the pre- 
vious results must be related to more subtle correlation 
of higher order nature. 

Lastly we would like to discuss the importance of early 
stages of the growth. (F^) might be written in the fol- 
lowing way 



(^r ) > = <n( 1 + ^) Q ) 



(51) 



(see Eq. (pi|)). Neglecting the correlations in time in the 
above product one may approximate 



di( i +w Ri n« i +^) a ) 



Numerical evaluation of the two objects in Ec 
that they differ by a few percent (see Fig. | 
merics indicate the scaling laws 



(52) 



( |52j ) shows 
. The nu- 



the initial and late stages of the growth are much more 
important than local correlations in the late stages. 

We checked also two-point time correlations (A„A„_fc) 
for some fc's. The results are plotted in Fig. ^. As it 
turns out from this graph, (A n A„_fc) ~ (A„) (A„_fc) up to 
statistical fluctuations. 



V. SUMMARY AND DISCUSSION 

The language proposed by Hastings and Levitov ap- 
pears to offer many appealing features. It generates DLA 
clusters in such a way that the conformal map from 
the circle to the boundary of the cluster is known at every 
instant. In this paper we examined carefully the numeri- 
cal procedure used to generate the conformal maps, and 
pointed out the advantages and the shortcoming of the 
algorithm. 

The new results of this paper pertain to the scaling 
behavior of the Laurent coefficients |Ffe| of the conformal 
map <&( n * ) and of the moments of A n which are related 
to moments of the field. We presented a theoretical dis- 
cussion of the exponents characterizing moments of \Fk\ 
and A ra . We pointed out the relations to the multfifractal 
analysis of the harmonic measure, and derived scaling 
relations. Of particular interest is the scaling relation 
Z?3 = D/2 that was derived first by Halsey and which 
appears here as a very natural consequence of the for- 
malism. 

One important result which is not adequately inter- 
preted in this paper is the sharpness of the distribution 
of F\ . This coefficient is proportional to the radius of the 
cluster, and its sharpness is directly related to the exis- 
tence of a universal fractal dimension independently of 
the details of the shape of the cluster. Understanding the 
sharpness appears to be connected to understanding the 
existence of universal fractal dimension, and we believe 
that it poses a very worthwhile and focussed question for 
the immediate future. 
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<n(l + Ai) a ) = cA 



o" 



2/D 



n<(! + A a ) =ciA 7 



2/D 



(53) 
(54) 



where ci/c > 1.06. 

To get further intuition we checked also the object 

k n 

(II( i+a o b > n ((i+w 

i=l i=k+l 

for various values of k. The results are shown in Fig. |^. 
As it seems from this graph, time correlations between 
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APPENDIX A: CONSEQUENCES OF THE 
ONE-FOURTH THEOREM 

In this appendix we prove that every univalent func- 
tion of the type ( pT| ) is bounded in a circle of radius 4Fi . 
This fact is based on two basic properties of univalent 
functions lOl: 
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1 . There is one-to-one correspondence between univa- 
lent functions of the form f(w) = a±w + a-iw + 
■ ■ • (5-class) and univalent functions of the form 
g(w) = aiw + a-i/w + a-2/w 2 + ■ ■ ■ (S-class). This 
correspondence is given by 



j(w) <-> g(l/w)~ 



(Al) 



2. The Koebe One-Quarter Theorem. The image of 
the unit disc under every function of class S con- 
tains the disc z : \z\ < 4|ai|. 



Consider a function of the form (11). This is a E- 

class function with l inea r coefficient F\ . Let us denote its 
conjugate (By Eq. (Al)) 5-class function as P(w). The 
linear coefficient of P is 1/F\. Consider now the small- 
est circle in the z— plane which bounds the ima ge o f the 

\z\ = R}. From Eq. [Al] it is 
\z\ = l/R} is the largest circle 



unit circle under {z 
clear that the circle {z : 
which is contained in the image of the unit disc under 
P. Thus the Koebe One-Quarter Theorem ensures that 
1/-R > 1/(4^1), which implies R < 4F X . 



APPENDIX B: ESTIMATE OF THE SCALING 



BEHAVIOR OF {|F^| 2 > 



To estimate the large k and large n dependence of 
F^jj , the components are first written as integrals over 
the boundary of the cluster 

F}$ = (l/2vr) / & n ~)(e ie )e ike dO 
Jo 



where 



z(s)e lke ^\E(s)\ds 



0(s) = / \E(s')\ds'. 
Jo 



(Bl) 



(B2) 



For the purposes of Sect. 4 we are interested in \Fl 



(™)|2. 



IF, 



(»»)|2 



L r L 



JO 



z(s)z*(s')e ik ^- e ^\E(s)\\E(s')\dsds' 

(B3) 

For a given value of n (or equivalently, of R ~ n 1 ^ ), 
an examination of Eq (B3) shows that for large k the 
fluctuations in the values of the integrands depend more 
crucially on the phase variations than on the field and 
radius variations. The phase varies appreciably when 
changes an amount 



AO w (7r/2Jfc). 



(B4) 



and therefore it is useful to split up the integral Eq.(Bl) 
into a sum of essentially independent contributions com- 
ing from the electric field singularities with exponents a. 



This exponent is determined by the scaling law relating 
the measure (which is proportional to AO) of a box to 
its size (As) a : AO - ((As) a /R) a [§. The integral is 
split into contributions made of contour sections of dif- 
ferent lengths (As) a dependent on the singularity but 
each giving rise to the same change AO. If one can esti- 
mate both the magnitude of the contribution of a specific 
multifractal electric field singularity a to the integral and 
the number of such contributions Af a (k,n) then one 
can write 



<|FM| 2 )~£|/ a ,(M)| 2 



J daAf a (k,n)\I a (k,n)\' 



where 



I a (k,n) = / \E(s)\z(s)ds 

J(AS) a 



(B5) 



(B6) 



To estimate Af a (k, n) we recall that by definition 

r s+As 

AO = I \E(s')\ds' - (As/R) a . (B7) 



^Frorn Eq.(|BJ) 



(As) a - i?( 7 r/2fc) 1 / Q . 



(B8) 



Using Eq.(B8) and the fact that the fractal dimension of 
singularities of size a is f(a), we can now also estimate 
the number of singularities of size a which contribute to 
the integral as 



Afjk) ~ (R/(As) a ) f ^ ~ (2fc/7r) /(Q)/tt 



(B9) 



To estimate \Ia(k, n)\ 2 we note that the major contribu- 
tion to Eq. (|B6| ) comes from the support of the harmonic 
measure where \z(s)\ w R. Accordingly 



\Ia\ 



R 2 



\E(s)\ds 



(As) c 



R 2 (A0) 2 - R 2 (2n/k) 2 , 
(BIO) 



where we have made use of Eqs.(B2), (B4). Combining 
Eq.(B5) with the estimates (B9) and (BIO) then yields 



(IF^I 2 ) ~ (R/Ak) 2 / da(2k/ir)^/ a 



(Bll) 



We note that the approximation adopted in this appendix 
differs from the delta-function assumption ( |35| ) in assert- 
ing that for high values of k the variation of the phase 
dominates the decay of the integrand compared to the 
rapid decorrelation of the field. One would guess that for 
k of the order of unity the field decorrelates faster due to 
the rapid variation over the arc length. For high values 
of k the phase decorrelation is strongly amplified and we 
adopt the assumption used here. 
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FIG. 3. The scaling exponents of the Laurent coefficients: (|F_ fc | 2 ) ~ n Xk . The values are obtained by averaging 400 
independent realizations of 10,000 particle clusters. 




FIG. 4. The scaling of (|^o n '| 2 ) (thick lines) and the sum of diagonal terms (thin lines, see Eq.(51)) with size n. Clearly the 
two have different scaling exponents. The solid lines are averages over 400 clusters of size 10,000, the dashed lines are averages 
over 30 clusters of size 100,000. 
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n 



FIG. 5. Scaling of the moments (A* ) with powers of n. The curves from top to bottom correspond to q — 0.5, 1, 1.5, 2, 2.5, 3 
and 3.5. The exponents —2qD2 q +i/D are in agreement with theoretical predictions (see text) and with numerical values for 
the generalized dimensions in the literature. 




n 



FIG. 6. The rescaled standard deviation aj™ of the Laurent coefficients of the map (see definition in text). For k ^ 1, 
cr^ n ' fluctuates around unity, corresponding to broad distributions. For k = 1 it tends to zero as n — ► oo, demonstrating the 
asymptotic sharpness of the distribution of F\. The solid lines are averages over 400 clusters of size 10,000, the dashed lines 
are averages over 30 clusters of size 100,000. 
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FIG. 7. Time-angle correlations of the field. In order to reduce statistical noise, the values plotted are averaged in bins 
[n, I. In]. 



1.06 




FIG. 8. The ratio of Fi approximated by neglecting time correlations and the full Fi: J\" =1 {(1 + ^i) a ) / (YYi^ii^ + Ai ) a ) 
(thick line). The quantities (fl^LiC 1 + A 0°) H7=k+A( 1 + A 0°> arc also plotted for k = 10, 100 and 1000. 
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FIG. 9. Correlations of the field. In order to reduce statistical noise, the values plotted are averaged in bins [n, l.Oln]. 
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